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The propagation of weakly nonlinear dust-acoustic waves in a dusty plasma containing two ion species 
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O ' with different temperatures is explored. The nonlinear equations describing both the quadratic and 

o : 

CN ■ cubic plasma nonlinearities are derived. It is shown that the properties of dust-acoustic waves depend 

O '. 

4) , substantially on the grain size distribution. In particular, for solitary dust-acoustic waves with a 

o 

\ positive potential to exist in a plasma with distributed grain size, it is necessary that the difference 
between the temperatures of two ion species be large that that in the case of unusized grains. 

1 INTRODUCTION 

Dust-acoustic waves are a phenomenon inherent in plasmas containing, besides electrons and 
ions, heavy (about 10~ 14 g) dust grains, whose charges are determined by the balance between the 
ion and electron currents onto on the grain surface. The processes occurring in dusty plasmas have 
been actively studied over the past two decades. The presence of dust in plasma not only modifies 
the conventional plasma modes, but also gives rise to new types of low-frequency electrostatic plasma 
oscillations, such as dust-acoustic waves [Tj, dust Coulomb waves [2], and dust-lattice waves [3]. A 
' considerable number of studies have been devoted to investigating linear and nonlinear properties 

>< : 

. of these waves. Weakly nonlinear small-amplitude waves are usually analyzed by expanding the 
equations in power series in a small parameter and using the method of coordinate stretching in 
accordance with the form of the dispersion relation [4], while large-amplitude waves are studied 
using the quasi-potential method [5] or gasdynamic approach |6J. 

Dust-acoustic waves propagate with a very low phase velocity (from several tens centimeters 
per second up to a few meters per second). Therefore, it can be assumed that the electrons and 
ions, whose thermal velocities are much higher than the wave phase velocity, are distributed in 
the wave field according to the Boltzmann law. A number of papers were devoted to studying the 
properties of nonlinear dust-acoustic waves in a plasma containing two ion species with different 
temperatures [7]— [1Q|. Such a situation can occur when two plasmas with different temperatures are 
mixed with one another. The time during which ions relax to a state characterized by one temperature 
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is usually much longer that the period of dust-acoustic oscillations (see |8J). Therefore, such a two- 
temperature plasma can be described by Boltzmann distributions for each ion species. It is also 
assumed that the electron temperatures of two plasmas are equalized very rapidly by virtue of the 
condition r ee /r n pa (m e /m i ) 1 / 2 (T i /T e ) 3 / 2 <C 1, where r ee and r n are the relaxation times of electrons 
and ions to equilibrium states; m e and rrii are the masses of an electron and ion, respectively; T e and 
Tj are the electron and ion temperatures. 

In [7], the propagation of a nonlinear dust-acoustic wave in a dusty plasma with two-temperature 
ions and a constant grain charge was considered. It was shown that both compression and rarefaction 
solitary waves, as well as double layers (waves having one nonzero asymptotics at infinity and 
described by the hyperbolic tangent function), can propagate in such a medium. In [8], it was 
noticed that it is necessary to take into account variations of the grain charge in the wave field, 
because the contribution of this effect to the nonlinear coefficient in the Korteweg-de Vries (KdV) 
equation, describing the propagation of weakly nonlinear dust-acoustic waves, is comparable with 
the hydrodynamic nonlinearity. 

In previous studies on the propagation of dust-acoustic waves in a two-temperature plasma, it was 
assumed that that dust grains are equal in size. However, it is well known that both laboratory and 
astrophysical dusty plasmas contain dust grains of different size [H]-[l4]. The distribution function 
of dust grains over their radii a can be described by the power law 

f(a) = Ga~\ (1) 

where G — is the normalizing factor. In calculations, it is usually assumed that the value of the 
exponent v is in the range 3.5 - 4 [15]; however, one can find in the literature alternative numbers. 
Thus, in [H] , the value of v was assumed to be in the range 0.9 - 4.5, and, in [17] , it was reported 
that, for Saturn's G rings, this parameter can even reach a value of 6-7. The model implies that the 
grain radii lie within the interval between the minimum value a,\ and the maximum value 02; it is 
assumed that, outside this interval, f{d) = 0. 

In the present paper, the propagation of dust-acoustic waves in a two-temperature plasma is 
studied with allowance for both adiabatic variations in the grain charge (caused by variations in the 
densities of the ambient electrons and ions in the wave field and, accordingly, variations in the fluxes 
of negative and positive charges onto the grain surface) and changes in distribution function (1). The 
size distribution function alters both the linear wave phase velocity (as compared to the case of a 
dust with unisized grains) and the nonlinear coefficient in the KdV equation. An important point is 
that this coefficient can change its sign. Previous studies have shown that it can take both negative 
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and positive values, which corresponds to compression and rarefaction dust waves, respectively. The 
latter is possible when the ratio between the densities of cold and hot ions is moderate, whereas their 
temperatures differ significantly. The goal of the present study was to investigate the propagation 
of solitary dust-acoustic waves in a plasma with two-temperature ions under the assumption that 
dust grain sizes obey distribution (1), to find the conditions required for the existence of waves with 
a negative and positive potential, and to compare the results obtained with results of calculations 
performed for a dust with equal-size grains. 

The paper is organized as follows. In Section 2, a set of basic equations is derived and dimensionless 
variables are introduced. In Section 3, the KdV equation and modified KdV equation (mKdV) with 
coefficients that take into account the distribution over dust grain sizes are obtained by expanding the 
basic equations in power series in a small parameter and using the method of coordinate stretching. 
In Section 4, the nonlinear coefficients of the equations obtained are analyzed. The main results are 
summarized in Section 5. 



2 BASIC EQUATIONS 

Let us consider a dusty plasma consisting of electrons, ions, and heavy dust grains of different 
size and mass. We assume that the dust contains N different species with the densities n j0 , charges 
Zjo, and masses rrij. Then, the condition of plasma quasi-neutrality can be written as 

N 

n ih0 + n m = n e0 + ^2 z jonjo, (2) 

3=1 

where n^ch^o, and n e o are the unperturbed densities of high- and low-temperature ions and elect- 
rons, respectively. We assume that the wave phase is much larger than the dust thermal velocity and 
ignore the influence of the thermal pressure on the dynamic of the dust component. In this case, the 
continuity equation and equation of motion for each dust species have the form 

+ Tx (n,u,) =0, (3) 

dv,j diij Z~e dd> , . 

at ox rrij ox 

where rij and Uj — are the density and velocity of the jth dust component and Zj — is the charge of 
its grains. This set of equations should be supplemented with Poisson's equation 

d 2 d) N 

— = 4vre(^ ZjUj + n e - n ih - n u ), (5) 

3=1 



in which the densities of electrons and both ion components obey Boltzmann distributions, 
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n ih = n iho exp (^jr^j > ( 7 ) 

nil = n uo exp (^-f^J ■ ( 8 ) 

Here, T e , T ifo , and T a are the temperatures of electrons and high- and low-temperature ions, respective- 
l.v. 

Variations in the charge of dust grains of the jth species are described by the balance equation 
of currents onto the grain surface, 

e Ht e + ih + ih ^ 

where 

I e = — eTTa 2 A%T e / m e ) l l 2 n e exp 



I ih = eKa 2 (%T ih /mi) 1/2 n ih I 1 



T 



ih 



eT:a 2 Ami/m i ) 1 ' 2 n il ( 1 - ^ 



$, = —Zje/cij is the surface potential of grains of the jth dust component and aj is their radius. 
The charging frequency of dust is determined by the expression u ch = ajio^j '(\/27rAi>j) , where uj pi = 
(4:%e 2 mio/mi) 1 / 2 is the ion plasma frequency and Xm — (7i/47re 2 n i o) 1 ^ 2 is the ion Debye length. The 
time during which a wave passes through a given point of the medium is = \/v p h ~ (A/ Xm)uJp^ , 
where u pc i = {A r KZ 2 e 2 ndo/m c i) 1 l 2 and A is the wavelength. We assume that the charge of a dust grain 
in the wave field varies quasi-statically, so the current balance equation has the form 

I e + I lh + In w 0. (10) 

Obviously, for this to occur, it is necessary that v ch T d ^> 1. Assuming that aj ~ 10~ 4 cm, Tj ~ 1 eV, 
Z ~ 10 3 , md ~ 10~ 14 r, and Zrido ~ n io, we obtain the following condition under which the dust 
grain charge in a long- wavelength dust-acoustic wave A^ A can be described adiabatically [T8] : 
n i0 3> 10 6 (X Di /X) cm -3 . For this condition to be satisfied, the plasma should be sufficiently dense. 
Otherwise, it is necessary to take into account the nonadiabatic character of grain changing in the 
wave field, due to which the initial perturbation transforms into a dust-acoustic shock wave [T9] . 
It follows from the above estimate that the model in which the grain charge is assumed to vary 
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adiabatically mainly applies to laboratory plasmas, while the use of this model to describe dust- 
acoustic waves in astrophysical plasmas, e.g., in Saturn's F rings ([ID]), is incorrect. 

In equation of dust motion (4), we ignored the frictional force. In dusty plasmas, the momentum 
loss of dust grains due to collisions with ions is usually negligibly small and can be ignored. Actually, 
the applicability condition of this approach is x = miTiiofid/ ' (rn^rido^pd) 1, where Uid = 
A7TAZ 2 e 4 ndo/my 2 T^ 2 is the collision frequency of ions with dust grains and A is the Coulomb 
logarithm. For the above parameter values and n^o ~ 10 6 cm -3 , we have x ~ 1CT 4 . In laboratory 
plasmas, the density of neutral particles is, as a rule, several orders of magnitude higher than the 
plasma density. In our problem, the condition under which the friction of dust by the neutral 
component can be ignored has the form fdn/^pd = 8\^a 2 m n n n VTn/3rridUJpd <C 1 [2Q], where n n 
is the density of neutral particles, m n is the mass of a neutral atom, and vm is their thermal velocity. 
For n d Q, we find that the density of neutral particles should be n n <C 10 15 cm -3 , which is usually 
satisfied with a large margin. 

Let us introduce the total dust density N = Ylf=i n jo an d the average charge Z = 
= J2f=i n joZjo/N, average radius a = 5^j=i a j n jo/N , an d average mass m = Y^j=i m j n jo/N of a 
dust grain. We also introduce the quantities 

ZoNT e T ih Tn 



T, 



eff 



neoTihTu + nihoT e Tu + nu T e T ih 
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^Dd = (T eff /ATcZ e 2 N)^ 2 , u pd = (AttZ e 2 N/m) 1 / 2 , and c D = (Z T eff /m) 1/2 , and normalize the dust 
density rij to N; the dust grain charge Zj to Z ; the grain mass rrij to m; the electron density n e0 
and ion densities riuo, to Z N = Ylf=i n joZjo'> the time t to 0J~ d x , the spatial coordinate x to Xod] 
the velocity Uj to cp; and the electrostatic potential and surface potential $j to T e ff/e. Then, the 
initial set of equations takes the form 

^ + £(»«)=0, (ID 

(12) 

at ox rrij ox 



d 2 cj) N | 

= Z i n i + s +6 _ 1 ( ex P(/ j i g< ^ ) ) - J i exp(-s0) - 5 2 exp(-/3s0)), (13) 
j=i 1 2 

expf/^s^ + $j)] = «i5iexp(-s0)(l - s$j) + a 2 5 2 exp(-/3s0)(l - f3s$j), (14) 

where the following notation is introduced: 5\ = nuo/n e o, 5 2 = niho/n e o, (3\ = Tu/T e , (3 2 = Tih/T e , 
(3 = T u /T ih = fa/ fa, a x = (A/^) 1 / 2 , a 2 = ((3 2 /fi) 1/2 , \i = rm/m e « 1836, and s = T eff /T u = 
(5 1 + 5 2 -l)/(fa + 5i + p8 2 ). 

In deriving Eq. (13), we used the condition Ylj=i n joZjO = 1) written in dimensionless variables. 

5 



3 WEAKLY NONLINEAR APPROXIMATION FOR DUST- 



ACOUSTIC WAVES 

We will search for the solution to the set of Eqs. (11)— (14) in the form 

Y = Y + eYi + e 2 Y 2 + e 3 Y 3 + . . . , (15) 

where Y = (<f),nj,Uj,&j) and Y = (0, nj , 0, $jo), by using the conventional method of coordinate 
stretching: 

^ = £ ^( X -V t), r = e 3/2 t. (16) 

In the zeroth order in e, we obtain the equation determining the surface potential of dust grains 
in an unperturbed plasma, 

exp(/?is$ i0 ) = ai£i(l - s$jo) + a 2 6 2 (l - Ps$ j0 ). (17) 

Since, in this equation, the surface potential $jo of different species of dust grains does not depend 
on j, the subscript j in <E>j will be further omitted. 

In the first order in e, the set of equations has the form 

r dun dun , . 

T dun Z~n dd>i , . 

AT TV 

Zj-iWj-o + ^2 Zjoriji + 0i = 0. (20) 

Integrating Eqs. (18)-(20) with the boundary conditions 0i, riji, Uji, Zji — > at £ — > ±00, we 
obtain the following relationships for the lowest order quantities: 



n n = -2-0i, Uji = —0i, Z,-i = 7i0i, (21) 

rrijVo rrijVo 

where the factor 71 is defined in the Appendix. For the linear wave velocity, we have 



In the second order in e, the set of equations takes the form 



- + + n j0 -^ + -^(n jlUjl ) = 0, (23) 
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du j2 duji duji Z j0 902 Z jl d<f>i 

V °~xF~ + ~a r ""i 1 ^ at ~ U ' 

9£ or o£ rrij o£, rrij o£ 



d 2 (j) 



r = 2^ Z i 2n io + 2^ + 2^ ^' on ^ 2 + 2 (5 1 + <5 2 - 1) ^ 



j'=i 



(24) 



(25) 



^2 = 7102 + 720?, (26) 

where the factor 72 is defined in the Appendix. 

Expressing drij 2 /d^ from Eqs. (23) and (24) and using Eq. (21), we obtain the KdV equation for 
a weakly nonlinear dust-acoustic wave, 



1 - + B^ + A?^ = 0. 



The coefficients A and B in this equation are 



(27) 



Vn 



2(1 + 7i) 



, 5M = 



s\5 l +(3 2 5 2 -(3 2 1 ) Zf n j0 3ti A ^0^0 

V n 4 ^ m 2 K 2 ^ 



<&1 + <&2 - 1 



'0 7=1 "*j 



i=1 j 



272- 



To analyze the case in which the ratio between the coefficients in Eq. (27) is zero, B/A = 0, we 
use the method of coordinate stretching, 



£ = e(x-V t), 



r = e 6 t. 



(28) 



Then, the time derivative drops out from Eqs. (23) and (24) and we obtain the following expressions 
for the second-order quantities: 



n j2 



u j2 



mj V 2 



rnjVo 



-Z j0 (j)2 



-Zjo4>2 



1 fzz% 



2 \ mj V 2 
W z h 



7i <Pi 



2 KrrijVg 



- 7i <Pi 



In the third order in e, the set of Eqs. (1 1) (14) takes the form 

dn j3 dn n du j3 d 

~ Vo ~dt + + nj0 ~dt ljc njlUj2 + nj2Ujl > = ' 

du j3 dUjt d Z j0 903 Zji 902 Z j2 901 



Vr 



d£ dr d£ 



rrij <9£ rrij <9£ rrij <9£ 



0, 



(29) 

(30) 
(31) 



9 2 0i 



s 2 ((3 2 - 5 1 -(3 2 5 2 ) 



j'=i 



j=i j=i j=i 

+ ^7^ 1 ^ 2 + ^3, 

0\ + o 2 — 1 

Zj3 = 7103 + 2720102 + 730?, 



(32) 
(33) 
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where the factor 73 is defined in the Appendix. 

Excluding dujs/d^ from Eqs. (30) and (31), differentiating expression (32) over £, and using Eqs. 
(21) and (29) and the condition B/A = 0, we arrive at the mKdV equation, 



where 



CIA 5 \- ^jo^jo 37i 4,-qWjo 7 1 n j0 
N - ,3^3 



4 72 Z j0 n j0 s 3 ffi + ^1 + /3 3 (5 2 



Under the condition 0i (77 = £ — -ur — > ±00) = , Eqs. (27) and (34) have the well-known solutions 
in the form of solitary waves. Their properties, such as the amplitude, the propagation velocity, and 
the sign of the electric potential, are determined by the values and signs of the coefficients entering 
into the nonlinear equations. Below, we will analyze the dependence of these coefficients on the 
plasma parameters. 



4 ANALYSIS OF NONLINEAR COEFFICIENTS 

Taking into account distribution function (1), we replace the sum over dust species with the 
integral, 

N ra 2 

J2n j0 = G a~ v da = N. (35) 

Then, averaging all the quantities under the summation symbol in the expressions for Vq, B/A and 
C/A, we obtain 

d 1 +d 2 -l Vq Vq 

n/A _5a 3 7 ic 2 7l 2 c 3 4 72 c 4 s 3 Pf+Jh+JPh 

C/A -2Vj~^f + 2V! + ^vJ~ l3 ~J 5 1 + 5 2 -l ' (38) 
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where 
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c 2 
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a 2 
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2-v n 2-v\2 
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(9 — 1A 3 
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n- 2 -A 
a 2 ) 


(4 - 


— is ) yz -|- v ) 




[a 2 — a x )° 
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V) 


\ a 2 


a l ) l a l 


a 2 ) 


f4 




(n 2 ~ 
\ a 2 


a l )\ a 2 


a l ) 






\ a 2 


a l ) l a l 


a 2 ) 


f4 


7/WzL -i- iA 




(a 2 - a x ) 




(2- 


-v)\\-v) 


(of 


- atrial - 


- a 2 3 -n 


(4- 


-u) 2 (3 + u) 


(al- 


»-a 2 mal-»- 






(1-u) 2 




- a 4 r)(al 2 - - 


a-2 2 -n 


(4 


-u)(2 + u) 




- a\~n 2 


1 



c 3 

I 4 - I" V) 

c 4 = b 2 . 

In deriving the above expressions, we took into account that the density rijo, charge Zjq, and mass 
rrij are dimensionless quantities; therefore, in taking integrals, we passed to dimensional variables in 
accordance with n — > n/N, Z(a) — > Z(a)/Z Q , and m(a) — > m(a)/m. We also used the dependences 
Z(a) = ($o/e)a, and m(a) = (47rp/3)a 3 to find expressions for the average values of Z and m: 

Zn = : 1 — , m = : ; — , 40 

e 2 - v a\- v - a\- v 3 4 - v a x 2 v - a x x ~ v 

using which expressions (39) can easily be derived. 

Formulas (39) allow a limiting transition to the case of a dust with unisized grains at a 2 — > a±, or 
v — > oo. In the former case, the uncertainty is eliminated by setting a 2 = a\ + Aa with Aa <C a and 
keeping the summands of the first order of smallness. In the latter case, the contribution from the 
distribution function is significant only in the vicinity of a± , irrespective of the width of the interval 
[di — a 2 ]. In both cases, coefficients (39) are reduced to unity and nonlinear coefficients (37) and (38) 
transform into the expressions obtained in |8j. 

The results of numerical analysis of coefficients (37) and (38) are presented in Figs. 1-6. Figure 
1 shows the dependence of the coefficient B/A on 5 2 for two values of Si for the case of a dust of 
unusized grains. This coefficient can vanish when the density of cold ions is smaller than the density 
of hot ions and the temperature ratio Tih/Tu is high enough. As the absolute value of the parameter 
5i (which characterizes the density of cold ions) increases, the zeros of the function B/A{5\) shift to 
the right. Figure 2 illustrates the dependence of the coefficient B/A on S 2 for different ratios a 2 /ai 
and a fixed value of the exponent v. It is seen that, as the width of the interval of dust grain radii 
increases, the critical point B/A = disappears. Figure 3 presents a series of plots B/A(8\) a 2 /ai and 
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different values of the exponents v. As was mentioned above, the case of large v values corresponds to 
a dust with grains nearly equal in size. For the given parameters, the coefficient B/A at v — > oo can 
be either negative or positive. However, as the exponent v decreases, the influence of the distributed 
grain size increases, the critical point disappears, and the coefficient becomes negative defined. Figure 

4 shows the dependence B/AiS-i) for different ratios between the temperatures of cold and hot ions. 
It is seen that the coefficient is very sensitive to the ion temperature ratio. For the given (wide) 
interval of dust grain radii, solitary rarefaction waves can exist only when the temperature ratio is 
larger than 35. At the same time, in a dusty plasma with equal-size grains, rarefaction waves can 
exist when the temperature ratio is 20, the other parameters being the same (see Fig. 1). 

The curves in the plane (Si — 5%) in Fig. 5 correspond to the points at which the coefficient 
B/A vanishes. The lower and upper curves correspond to the left and right zeros of the function 
B/A. For 5\ values in the range [0 — 0.85], there are two points on curve 1 at which the nonlinear 
coefficient vanishes. The existence of two critical points of this coefficient in a dusty plasma with 
two-temperature ions was revealed in [7]. At Si = 0.85, the critical points merge, and, at Si > 0.85, 
the coefficient is negative defined. The region in which B/A > (rarefaction dust wave) lies inside 
the sector bounded by the curve. The larger the spread in dust grain sizes, the narrower the domain 
of existence of rarefaction waves (see curve 2). At v < 2.5 (the other parameters being the same), 
the coefficient B/A is negative defined and nonlinear dust-acoustic waves are compression waves. 
The dependence C/A(Si) under the condition B/A = has the shape of a loop (see Fig. 6). The 
two points at which the coefficient B/A vanishes correspond to two different curves C/A. Curves 1 
and 2 refer to the left and right zeros of B/A, respectively. At Si ~ 0.85 , the zeros of the coefficient 
B/A merge; accordingly, the values of C/A also merge at that point. It should be noted that, on 
curve 2, the coefficient C/A can take negative values. This means that double layer structures can 
form in a dusty plasma. We also note that, at C/A = (curve 2) , the second-order nonlinearities 
are mutually cancelled. In this case, the nonlinear properties of a solitary wave are determined by 
the third-order terms and can be analyzed in a similar way [7]. 

5 CONCLUSIONS 

As the range [01,02] of dust grain sizes extends and/or the exponent v of distribution functions 
(1) decreases, the positive values of the nonlinear coefficient B/A n the KdV equation decrease and 
its negative values increase in magnitude (see Figs. 2, 3). In this case, the amplitude of a rarefaction 
solitary wave increases and that of a compression solitary waves decreases as compared to those in a 
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dusty plasma with unisized grains. Even if the ratio between the temperatures of cold and hot ions, 
Tu/Tih, , is quite sufficient for the development of solitary waves with a positive potential in a dusty 
plasma with unisized grains, it can be insufficient for their development in the case where the grain 
sizes obey distribution (1) (cf. curve 1 in Fig. 1 with curve 3 in Fig. 4). Thus, in order for rarefaction 
solitary dust-acoustic waves to exist in a plasma with two-temperature ions and distributed dust 
grain size, it is necessary that the temperature ratio between two ion components be larger than that 
in a dusty plasma with unisized grains. 
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the Russian Federal Agency on Education (project no. RNP 2.1.1.3483) and the Russian Foundation 
for Basic Research (project nos. 05-02-17070 and 06-02-16819). 



Let us define constants relating perturbations of the dust grain charge to the wave potential 0. 
We expand the equation 
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APPENDIX 



exp[/5is(0 + $)] = ai<5i exp(-s0)(l - s$) + a 2 8 2 exp(- / 9s0)(l - /3s<E») 



in power series in <f> and In the first order in order e, we have 



0i [ft exp[(3 lS %) + Mi(l - s$ ) + a 2 S 2 (3(l 



+$i[/3i exp[/3is$ ) + ai<Ji + a 2 5 2 p] 



0. 



Introducing the quantities 



fi = pi exp[/3is$o) + ai<$i(l - s$o) + a 2 5 2 /3(l - /3s§ ) 



r 2 = Pi exp[/3is<l>o) + cti^i + a 2 S 2 p, 



we obtain 



$i = -=^i. 



1 2 

Using the relationship Z\ = $i/<l>o for dimensionless quantities, we find 



r 



i 



Zi = 7i0i, where 71 = 



r 2 $ 
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The factor 71 is the same for any jth dust species; therefore, the subscript j in Qjk and Zj k (k = 1, 2, 3) 
can be omitted. 

In the second order in e, we have 

i\0 2 + r 2 $ 2 = -^ 1 2 sexp(/3 lS $ o )(0i + $i) 2 + «iM0i$i + (1 - s$ o )0?/2) + 

+a 2 5 2/ 3 2 s(0 1 $ 1 + (1 - /3s$ o )0 2 /2). 

Using the relationships $1 = 7i$o0i and Z 2 = $ 2 /$o, we obtain 

r 3 



Z2 = 7i02 + 720i, where 7 2 



r 2 <V 

r 3 = -^ 1 2 sexp(/3 1 s$ )(l + 7i$o) 2 + aiM 7 i$o + ( X ~ s$ o)/2) + 
+a 2 (5 2 /3 2 s(7 1 $ + (l- / 3s$o)/2). 

Finally, in the third order in e, we obtain 

Ti0 3 + r 2 $ 3 = 2r 3 0i0 2 + r 4 0?, 

^3 = 7103 + 27 2 010 2 + 730? 



where 



73 = T 4 = ~Pls{l + 7i$o) exp(/3 lS $ )(7 2 $ + (3 lS (l + 7 i$o) 2 /6)- 

1 2^0 

-« 1 (5 1 s 2 (7 1 $ /2 - 72$o/s + (1 - s$ )/6)- 
-a 2 5 2 /5 3 s 2 ( 7 i<f o/2 - 72$ o//3s + (1 - /3s$ )/6). 
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FIGURE CAPTIONS 

Fig. 1. Nonlinear coefficient B/A as a function of 82 for the case of a dust with equal-size grains at 
different temperatures of the cold ion component: (a) ft= 1) 0.03, 2) 0.05, 3) 0.07 (81 = 1, ft = 1); 
(b) ft= 1) 0.02, 2) 0.03, 3) 0.05 (<5i = 2, ft = 1). Key: 1. (a); 2. (b) 

Fig. 2. Nonlinear coefficient B/A as a function 82 for different values of the ratio 02/01: 1) 1, 2) 
1.5, 3) 2, and 4) 2.5 (81 = 1, ft = 0.04, ft = 1, 1/ = 3.5). 

Fig. 3. Nonlinear coefficient as a function of 8 2 for different values of v. 1) 7, 2) 6, 3) 5 
(8 1 = 1, ft = 0.04, ft = 1, a 2 /ai = 10). 

Fig. 4. Nonlinear coefficient £>/A as a function of <5 2 for different values of ft: 1) 0.02, 2) 0.025, 
3) 0.03 (S 1 = 1, ft = 1, a 2 /ai = 10, 1/ = 3.5). 

Fig. 5. Curves defined by the condition B/A = in the plane [Si — 82] for 1/= 1) 3.5, 2) 3 (ft = 0.03, 
ft = 1, a2 /ai = 10). 

Fig. 6. Coefficient C/A as a function of Si for £>/A = 0. Curves 1 and 2 correspond to left and 
right zeros of the coefficient B/A, respectively (ft = 0.03, ft = 1, 02/01 = 10, v — 3.5). 
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Fig.6 
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